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We investigate how the macroscopic spatial structure of broken inversion symmetry manifests in 
noncentrosymmetric superconductors, by the microscopic broken inversion symmetry of the crystal 
structure. Based on the time-dependent Ginzburg-Landau theory including the Pauli paramagnetic 
effect and the Rashba interaction, we demonstrate that the centrosymmetric structures of the inter- 
nal field and the screening current are broken macroscopically. The flow structure of paramagnetic 
supercurrent spontaneously induce the flux flow without applying external currents. 
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I. INTRODUCTION 

Recently, several heavy fermion superconductors 
breaking inversion symmetry in the crystal structure 
were successively discovered, such as CePtaSi,^ Ulr,^ and 
CeRhSia^ and attract much attension. In the noncen- 
trosymmetric superconductors, there appear some exotic 
properties through the spin-orbit coupling of the Rashba 
interaction which induces Fermi surface splitting between 
two helical spin projections.*'^ For the superconductivity 
on the Fermi surfaces, the special pairing form is only 
possible in the case of the spin-triplet pairing channel, 
or the upper critical field can exceed the Pauli-Clogston 
limiting field even for the spin-singlet pairing.^ Further- 
more, application of a magnetic field makes the super- 
conductivity more exotic, such as helically modulated 
superconducting order parameter,^ and spin polarization 
Ms oc n X Js induced by a supercurrent Js , where n is the 
unit vector along the polar axis of the noncentrosymme- 
try. Then, the paramagnetic supercurrent V x Mg also 
appears in association with MsA2ii£ 

The purpose of this paper is to show how the macro- 
scopic structures of broken inversion symmetry appear in 
noncentrosymmetric superconductors whose microscopic 
inversion symmetry is broken in the crystal structure, 
with a view to detecting experimentally direct evidence 
of the noncentrosymmetric superconductors. In the pres- 
ence of the time reversal symmetry at a zero field, we can 
not expect the broken-centrosymmetric spatial structure, 
since the band structure in momentum space still con- 
serves the inversion symmetry even in the noncentrosym- 
metric crystal structure. ^^ Under magnetic fields, by the 
combined two broken inversion symmetries of space and 
time, exotic noncentrosymmetric spatial structures will 
be revealed. We will see how these structures of the in- 
ternal field and the screening current appear both in the 
Meissner state and in the mixed state by the paramag- 
netic effect, using the simulation of the time-dependent 
Ginzburg-Landau (TDGL) theory. 

We describe the TDGL equation for the noncentrosym- 
metric superconductors and explain the method of our 
calculation in Sec. II. In Sec. Ill, we consider noncen- 



trosymmetric spatial structures of the screening current 
and the penetrating internal magnetic field near surfaces 
in the Meissner state. In Sec. IV, we study exotic spatial 
structures of the current and the internal field around 
a vortex. These structures induce the spontaneous flux 
flow without applying external currents. The last section 
is devoted to the summary. 



II. NUMERICAL SIMULATION FOR 

NONCENTROSYMMETRIC 

SUPERCONDUCTORS 

Numerical simulation of the TDGL equation is one of 
methods to investigate dynamical and static properties 
of superconductivity. The theory assumes that time evo- 
lutions of an order parameter A and a vector potential A 
depend on the relaxation processes of a superconducting 
free energy to the equilibrium. Thus, the TDGL equation 
is obtained bjii^ 
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in an appropriate unit for times, where i^ is a super- 
conducting free energy, and S/Sx indicates a functional 
derivative with respect to x. 

Hamiltonian of a noncentrosymmetric s-wave super- 
conductor with the Rashba interaction is expressed as 

7Y = dr lip^ i -—era + (ap x n + ^sB) • ct j -0 

-I- gipha-y'iljiphayip} , (2) 

with the electron momentum p — —iKV, the spin-orbit 
coupling constant a, the Bohr magneton /zb, the pair- 
ing interaction g, and the 2x2 unit matrix ctq, and the 
Pauli matrix a = (ca;, cr^, CTz). -0 is a field operator of 
the conduction electrons. For simplicity, we consider the 
case of a single component order parameter, neglecting 
the mixture of different parity components of the order 
parameteriSiii We also neglect the interband pairing be- 
tween split bands of ±p x n spin projections, since we 
consider the case when the band splitting energy due to 



the Rashba-interaction is enough large, as shown in the 
band structure calculation of CePt3Si:^i 

From the Hamiltonian ^, the GL free energy is ob- 
tained as^ii^ 
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for the order parameter ^, where a, b, and 7 are coef- 
ficients depending on the details of Fermi surface struc- 
ture, (/)o is the flux quantum. B and H are, respectively, 
internal and external magnetic fields. The last term with 
material dependent coefficient e (oc a^B) appears due to 
the absence of the inversion symmetry. The expression 
of Eq. © is also applicable to the single component case 
in p-wave pairing, as well as s-wave pairing. In a dimen- 
sionless form, Eq. Q is written as 



1 
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where A is the normalized order parameter, k is the GL 
parameter, Jg is the supercurrent density defined by 



Js = - (1 - T) Re [A* (-iV + A) A] /k^ 



(5) 



Here, e is normalized by ttooCo/^'/'o, where oq and ^0 are 
zero temperature values of a and coherence length, re- 
spectively. Using Eqs. 1^ and Q, we obtain the TDGL 
equation coupled with Maxwell equations 
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where the spin polarization Ms ex n x Jg is defined as 
1-T 



Ms = -e- 



-n X Re [A* {-iV + A) A] . 



(8) 



In the case of uniform B, we see from Eq. Q that 
A has phase factor of the helical state as A(r) — 
^(•j.)g-i£(nxB)r_7 -y^i^en this form of A is substituted to 
Eq. 10, the e-dependence in A* (— iV -|- A + en x B) A 
vanishes. Then, no additional diamagnetic current are 
induced by the phase factor of the helical state in the 
case of the uniform B. On the contrary, V x Mg in Eq. 
O gives the e-dependence due to the phase of the he- 
lical state, because en x B is absent in Eq. ||Sll. When 




FIG. 1: Spatial distribution of the current flow in the Meiss- 
ner state. We plot V x B in (a) and V x B — V x Bo in (b) 
for the region of 50^o x 50^o with open boundary. The length 
of arrow indicates the amplitude of the current. 



B has spatial variation, the phase factor of the helical 
state becomes nontrivial, and the e-dependence appears 
also from other terms in Eqs. (jHI-lIHll- Our calculations 
totally include contributions from the phase factor of the 
helical order parameter in Eq. ©, spin polarization M^, 
and effects of nonuniform B. These induce the noncen- 
trosymmetric structure in the macroscopic length scale. 
We perform numerical simulations following the time- 
evolution given by Eqs. ©-((SJ, as in previous works for 
centrosymmetric superconductorsiiS4iSiiLi& We consider 
the case n = (0, 0, 1) as in CePtaSi and CeRhSisfi^i^ and 
apply the magnetic field H = (iJ, 0,0) along the x axis. 
The calculations are done in a two dimensional square 
region parallel to yz plane. The region size is 50^o x 50^o 
(i.e., -25 < y/^o < 25 and -25 < z/^q < 25). Outside of 
the region we set A = and B = H. In our calculation. 



y=0 




FIG. 2; Profile of tlie internal magnetic field in the Meissner 
state along a plane parallel to n at y = 0. We plot B, Bo, 
B — Bo and Ms as a function of z. The inset focuses the 
behavior at the bulk region. 



near the surface at z = 25^o (^ = — 25^o), the screening 
current V x B is enhanced (suppressed) so that B is sup- 
pressed (enhanced) compared with Bq- This indicates 
that the field penetration is different between -|-n and 
— n-chrections.^^ 

In the closed system as in our case, due to the cur- 
rent conservation, the backfiow current to the param- 
agnetic current flows toward the opposite — y-direction 
far from the surface. Therefore, contraly to the region 
near the surface, the screening current V x B is sup- 
pressed (enhanced), so that B is enhanced (suppressed) 
compared with Bq in the inside region — 4^o < z < 19^o 
(-21^0 < z < -4^0). 

On the other hand, in the profile of B along a plane 
perpendicular to n at z = 0, B ^^ Bq, i.e., B has sym- 
metric profile for the inversion y <-> ~y. 



T — 0-5Tc^H=Q,e=o and k = 2.0, which is enough large k 
to study the qualitative properties in type-II supercon- 
ductors, while we use smaller k than that of CePtsSi for 
the reason of the system size and the efficiency of the 
simulation. We set e = 1 and qualitatively study the 
contribution of the e-dependent terms due to the non- 
centrosymmtry. 



III. SURFACE STRUCTURE IN THE 
MEISSNER STATE 

First, we discuss the screening current and the pene- 
trating internal field near surfaces in the Meissner state, 
when H ~ 0.04_ffc2,T=o,e=o- We show the flow of the 
screening current V x B in Fig. ^ (a), and the pro- 
file of the internal field B along a plane parallel to n 
at y = in Fig. El We see that, due to the screening 
current along the surface, B is dumped toward inside of 
the superconductor. We also calculate the internal field 
Bq when e = for the case of the conventional super- 
conductor with centrosymmetry. Then, to estimate the 
contribution of the e-dependent terms due to the noncen- 
trosymmetry in GL Eqs. (0)-®, we show the difference 
V x B - V X Bo in Fig. ^h) and B - Bq in Fig. H From 
this difference, we see that spatial structure of B breaks 
centrosymmetry in the direction of n. That is, B(r) is 
not symmetric for the inversion z <-> — z. Qualitatively, 
the noncentrosymmetric structure of B is due to the con- 
tribution of Ms appearing near surfaces at z = ±25^o- 
We also show Alg in Fig. |2] The contribution of e also ap- 
pears from other terms than Mg, such as the phase factor 
of the helical state, when B has spatial variation. Near 
the surface, Mg dominantly contributes to the sign of 
B — Bq. Far from the surface, there appear other contri- 
butions compensating Mg. As is shown in Fig. Q^b), the 
paramagnetic supercurrent V x Mg flows toward the same 
-|-y-direction near both surfaces at z — ±25^o- Therefore, 



IV. VORTEX STRUCTURE IN THE MIXED 
STATE 

Next, we discuss the vortex structure in the mixed 
state, when H = 0-15Hc2.t=o.s=o- In Fig- Ela): we show 
the flow of the current V x B, picking out 6^o x 6^o region 
around a vortex. While we set the position of the vortex 
center as (y, z) = (0, 0) in the figure, the winding cen- 
ter of the current does not coincide to the vortex center, 
shifting to the n-direction. This is due to the param- 
agnetic supercurrent V x Mg. As shown in Fig. Ofb), 
V X Mg flows toward the +y direction at the vortex core. 
And backflow of V x Mg returns to —y-direction outside 
of vortex. Figure^ shows the profile of the internal field 
and the current along the plane parallel to n through the 
nearest meshpoint to the vortex center. In Fig. Efa), we 
also show Mg and B — Mg in addition to B. The diamag- 
netic component B — Mg is almost symmetric for z <-!■ — z, 
as in the conventional centrosymmetric superconductors. 
The peak position oi B — Mg coincides to the vortex cen- 
ter. However, Mg shows almost antisymmetric structure 
for z ^^ ~z, since Jg of Mg oc n x Jg changes the flow 
direction at z ~ 0. Summing up S — Mg and Mg, we 
understand that B has noncentrosymmetric structure at 
the vortex. 

In Fig. EJb) , we plot the y-component of the current 
(V X B).y with (V x Ms)y and (V x B)y - (V x Mg)j,. The 
flow of the diamagnetic current component (V x B)j, — 
(V X Mg)y is winding around the vortex center. The 
slope of Mg in Fig. 2Ia) gives the paramagnetic current 
component (V x Mg)j,. Since Afg is increasing near z = 
in Fig. EJa), (V x Mg)j, flows +y direction at the vortex 
core. Therefore, the position where (VxB)j, = deviates 
from z = to -|-z-direction. The noncentrosymmetric 
vortex structure, breaking cylindrical symmetry, appears 
due to the antisymmetric spatial structure of the spin 
polarization Mg. 

On the other hand, in the profile along a plane perpen- 
dicular to n through the nearest meshpoint to the vortex 
center, Mg ^ and (V x Mg)^ r^ 0, so that B is sym- 
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FIG. 3; Spatial distribution of the current flow around a 
vortex. We plot V x B in (a) and V x Ms in (b) for a fo- 
cused region of 6^0 x 6^0 • The length of arrow indicates the 
amplitude of the current. The vortex center is located at 
iy,z) = {0,0). 
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FIG. 4: Profile of the internal magnetic field (a) and the 
current (b) around a vortex along the plane parallel to n 
through the nearest meshpoint to the vortex center. We plot 
B, Ms, B - Ms and |A| in (a) and (V x B)y, (V x Ms)y, 
(V x B)y - (V X Ms),; in (b). 
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metric and (V x B)^ is antisymmetric at z = for the 
inversion y *-^ —y, as in the centrosymmetic supercon- 
ductors. 

The flow structure of the paramagnetic supercurrent at 
the vortex gives similar effects as in the situation when an 
external current is applied toward +y direction. Then, 
in our simulation, the paramagnetic supercurrent spon- 
taneously induces the flux flow even in the absence of 
the external current. Figure |S1 shows the snapshot of the 
spatial structure of I A (r) I when H — 0.3iJc2,T=o,e=o;"^ In 
this simulation without applying external currents, new 
vortices are generated in the upper (positive z) surface 
continuously, and finally they move out to the lower (neg- 
ative z) surface after moving to lower side. The vortices 
do not form stationary lattice permanently. In our sim- 
ulation, the local voltage (ex; dA/dt)^ appears in the 
vortex. 

The above-mentioned broken-centrosymmetric fea- 
tures of the noncentrosymmetric superconductors ap- 
pear, when the magnetic field H is applied perpendic- 
ular to n. In the study of the vortex structure when H 
is parallel to nj-'^^ the spin polarization Ms appears to- 
ward the radial direction around a vortex, and the vortex 
structure remains to be cylindrical symmetric. 



V. SUMMARY 

We have shown how broken-centrosymmetric spatial 
structures appear in the macroscopic scale of the noncen- 
trosymmetric superconductors. We investigate the sur- 
face structure in the Meissner state and the vortex struc- 
ture in the mixed state, by numerical simulations of the 
TDGL equation. When an applied field H is parpendicu- 
lar to n, the broken-centrosymmetric spatial structures of 
the internal magnetic field and the screening current are 
induced by the characteristic paramagnetic effect due to 
the Rashba interaction. Namely, the field penetration in 
the Meissner state becomes anisotropic between +n and 
— n-directions. And the antisymmetric spin polarization 
Ms around a vortex breaks cylindrical symmetry of the 
vortex. The paramagnetic supercurrent at a vortex cen- 
ter given by the antisymmetric Ms may induce the flux 
flow without applying external currents. We hope that 
these phenomena will be observed in future experiments, 
as direct evidence of the noncentrosymmetric supercon- 
ductors. 



FIG. 5: Snapshot for the spatial structure of |A(r)| in the 
mixed state without applying external currentsi— 
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